import numpy as np
import matplotlib.pyplot as plt

# 定义 x 的范围
x = np.linspace(0.1, 10 * np.pi, 1000)  # 避免 x = 0

# 计算 y = sin(x)
y = np.sin(x)+10

# 过滤出 sin(x) > 0 的部分
positive_mask = y > 0
x_positive = x[positive_mask]
y_positive = y[positive_mask]

# 计算双对数坐标形式
log_y_positive = np.log(y_positive)
log_log_y_positive = np.log(log_y_positive)

# 绘制结果
plt.figure(figsize=(12, 6))

# 原始坐标系中的函数图像
plt.subplot(1, 2, 1)
plt.plot(x, y, label='y = sin(x)')
plt.title('y = sin(x)')
plt.xlabel('x')
plt.ylabel('sin(x)')
plt.legend()

# 双对数坐标系中的函数图像
plt.subplot(1, 2, 2)
plt.plot(x_positive, log_log_y_positive, label='Log-Log Plot of sin(x)')
plt.title('Log-Log Plot of y = sin(x)')
plt.xlabel('x')
plt.ylabel('log(log(sin(x)))')
plt.legend()

# 显示图像
plt.tight_layout()
plt.show()